Ultracold collision properties of metastable alkaline-earth atoms 
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Ultra-cold collisions of spin-polarized 24 Mg, 40 Ca, and 88 Sr in the metastable 3 P2 excited state are 
investigated. We calculate the long-range interaction potentials and estimate the scattering length 
and the collisional loss rate as a function of magnetic field. The estimates are based on molecular 
potentials between 3 P2 alkaline-earth atoms obtained from ab initio atomic and molecular structure 
calculations. The scattering lengths show resonance behavior due to the appearance of a molecular 
bound state in a purely long-range interaction potential and are positive for magnetic fields below 
50 mT. A loss-rate model shows that losses should be smallest near zero magnetic field and for fields 
slightly larger than the resonance field, where the scattering length is also positive. 

PACS numbers: 32.10.Dk, 33. 55. Be, 34.10.+X, 39.25.+k 



Exploring ultracold collision physics with metastable 
(nsnp) 3 P2 alkaline-earth atoms (ij [|, ||, |J, || promises 
insights that complement those obtained from the more 
conventional atomic species used in laser cooling. Unlike 
alkali-metal atoms |3], ^| the most-common alkaline-earth 
atoms have no nuclear spin, which greatly simplifies a 
theoretical description Furthermore, in contrast to 
metastable noble-gas collisions^ the low electronic en- 
ergy of the alkaline-earth metastable states ensures the 
absence of collisionally induced Penning or associative 
ionization. 

Atomic collisions in the ultracold regime play a crucial 
role on the road toward quantum degenerate gases. We 
show here that polarized metastable (nsnp) 3 P2 alkaline 
earth systems with projection m — 2 along the magnetic 
field might satisfy the key requirements for this quest: a 
positive scattering length and a favorable ratio of elas- 
tic to inelastic collision rates for certain magnetic field 
strengths. In addition, we show how a new type of pure 
long range molecular states can allow the resonant mag- 
netic field control of the scattering length. These states 
arise due to an interplay of an anisotropic (quadrupole) 
interaction and the magnetic field. Similar anisotropy 
has been predicted to exist for polar molecule collisions in 
an electric field jl^]. Many-body systems with anistropic 
interactions might now be explored [pd| . 

Trapping of metastable alkaline-earths in magneto- 
optical (|, H U> Hi and magnetic traps (|, [| [|, [l^] has been 
demonstrated. A ^1 s magnetic trap lifetime for densities 
near 10 10 atoms/cm 3 is observed with metastable 88 Sr in 
experiments at Rice and Tokyo university^, 13 . These 
trap lifetimes are limited by background collisions rather 
than by the radiative lifetime, which is on the order of 
minutes Q. 

Molecular potentials — We have calculated short-range 
adiabatic potentials correlating to two 3 P2 atoms using a 
molecular ab initio relativistic valence bond method, as 
previously described p3. A variety of attractive and re- 



pulsive potentials, and a number of short-range avoided 
crossings appear among them. Strong collisional loss pro- 
cesses are likely if the atoms approach one another 
on the attractive curves. The section below shows how 
a magnetic field can be used to manipulate long-range 
properties in such a way as to create pure long-range 
molecular states and prevent some of these losses. 

Atomic structure calculations and Raylcigh- 
Schrodinger perturbation theory allow us to express the 
long-range adiabatic potentials in terms of multipole 
moments of the atoms The matrix elements of the 
Hamiltonian for two 3 P2 atoms up to order 1/i? 6 , where 
R is the internuclear separation, are |lq| 
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Here \x) and \y) are product states | 3 P2, ^i)\ 3 P2, ^2), 
is the projection of the total electron angular momentum 
j a = 2 of atom a on the internuclear axis, and O = Qi + 
which is a good quantum number in the Hund's case 
(c) coupling scheme. The two terms on the right hand 
side of Eq. |l] describe the quadrupole-quadrupole (QQ) 
and the second-order dipole-dipole (DD) interaction. 

Matrix elements of the QQ interaction are expressed 
in terms of the quadrupole moment Q of the 3 P2 state [Q. 
For this paper the value of Q is recalculated using a 
more accurate many-body technique [0. The values 
are 8.46(8) E h <4, 12.9(4) E h a%, and 15.6(5) E h a% for 
metastable Mg, Ca, and Sr, respectively. Here Eh is 
a Hartree, 1 oq is 0.0529 nm, and one-standard devi- 
ation uncertainties, based on convergence studies of the 
many-body theory, are given in parenthesis. The Cq 3 are 
intermediate dipole-dipole dispersion coefficients, calcu- 
lated following Ref. |j, and tabulated in Table | The 
(x\Ajj'\y) only depend on angular momentum algebra 
and are given in Ref. Oil ], 

The Hund's case(c) adiabatic potentials that connect 
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TABLE I: Intermediate dispersion coefficients C" in units of 
10 3 EhdQ. The C% are symmetric in j and j'and have a 10% 
one-standard deviation uncertainty. 
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to our short-range valence-bond calculation are obtained 
by diagonalizing Eq. [11 within the | 3 P2, SI i } 1 3 i-2 , ^2) basis 
(See for examle Ref. jjl]].). 

Magnetic field — A magnetic field B modifies the long- 
range interaction of two metastable alkaline-earth atoms 
profoundly. For each atom a the approximate Zeeman 
Hamiltonian Hz a =A*s(ja+s Q )-B has been added to 
the other terms in the molecular Hamiltonian. Here, s a 
is the electron spin of atom a, /ib is the Bohr magneton, 
and we have assumed that the atomic gyromagnetic ratio 
of the electron orbital angular momentum and spin are 
one and two, respectively. The magnetic field lifts the 
degeneracy with respect to the projection m a of j Q along 
the space-fixed magnetic field direction. Forces due to the 
interatomic interaction break this space-fixed quantiza- 
tion, and align the molecular angular momentum along 
the internuclear axis. This results in loss of polarization 
of the angular momentum. 

The long-range adiabatic molecular potentials in the 
presence of a magnetic field are found by diagonalizing 
H Z i+H Z 2+Hir within the 3 P 2 + 3 P2 product basis. These 
potentials U (R, 9b) not only depend on the field strength 
but also on the angle 9b between the internuclear axis 
and B. Figure [l] shows all Sr2 gerade adiabatic potentials 
for a 10 mT field and 9b — 45°. The nine dissociation 
limits are separated by about E/Ub =10 mK and are 
labeled by M = mi + 777,2, where M = +4(— 4) for the 
highest (lowest) dissociation limit, ks is the Boltzmann 
constant. 

Figure [l] shows a multitude of avoided crossings for in- 
ternuclear separations where the Zeeman splitting is com- 
parable to the quadrupole-quadrupole interactions, that 
is, avoided crossings appear for R rj Rb = \J Q 2 /(^bB). 
For Sr* and B = 10 mT, R B is ~ 100 a . For R < R B 
the QQ interactions dominate and the Hund's case (c) 
Sl^ molecular potentials correlating to two 3 P 2 atoms are 
recovered. 

Experiments trap the "low-field seeking" spin- 
polarized m a = +2 state H 0]. The top-most potential 
in Fig. [j] dissociates to this limit. The radial and an- 
gular dependence of this potential t/ mi= 2,m 2 =2(P, #b) is 
shown in Fig. |^. Four long-range hills and valleys are vis- 
ible for R > Rb and a repulsive hard core that is nearly 
independent of angle appears for R m Rb ■ 

The behavior of U2, 2 (R,9b) is further illuminated 
when it is expanded in terms of Legendre polynomials 
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FIG. 1: Long-range gerade adiabatic potentials correlating to 
two 3 P 2 Sr atoms for B = 10 mT and 6 B = 45°. 
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FIG. 2: Polar plot of the adiabatic potential for two spin 
polarized Sr* atoms at _B=10 mT. Equipotential curves are 
shown as dashed, thick dashed, and solid contours for energies 
that are lower than, equal to, and higher than the dissociation 
energy of the two atoms, respectively. The contour values are 
given in mK. The arrow indicates the direction of B. 
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U2,2(R,9 b )= J2 U l {R)P l (x), 
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where x = cos9b- The Ul{R) of Sr2 are shown in Fig. ||. 
For R ^> Rb, U4 is largest and is dominated by the QQ 
interaction. In Fig. |^ it gives rise to the hills and val- 
leys. For R ps Rb, contributions of other components 
are appreciable, while for R < 60 an the repulsive Uq(R) 
dominates. The inset of Fig. || shows that Uq(R) has a 
shallow attractive region for large R, that becomes deeper 
and wider for increasing B. This so-called "Zeeman- van 
der Waals" well is a consequence of an interplay between 
Zeeman and anisotropic quadrupole-quadrupole interac- 
tions. 
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For R 3> Rb, an analytical expression for [72,2 
can be derived based on perturbation theory around 
| 3 P2, wi)| 3 P2, WI2) with the quantization axis along B. 
We find for the contribution of the QQ interaction 
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and for the DD interaction 
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The monopole coefficient Cg '' is independent of B and is 
1000 E h a%, 3300 £7^, and 6200 E h a% for Mg*, Ca*, and 
Sr* respectively. Consequently, the long-range behavior 
of Uo (R) is given by an attractive 1/R 6 potential that is 
independent of B. 
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FIG. 3: The coefficients Ul(R) of Eq. |2| for Sr2 as a function 
of internuclear separation and B =10 mT. The inset shows 
U (R) for B = 5 mT and 10 mT. 



Scattering lengths — A rigorous description of elas- 
tic collisions between spin polarized metastable 3 Pj — 
2,m = 2 atoms should start from B = theories! 19 
and extend those to include the Zeeman interaction^]. 
Such a model is beyond the scope of this paper. In- 
stead, we take advantage of our understanding of the 
adiabatic potentials in a magnetic field. Nuclear mo- 
tion couples the electronic potentials via so-called non- 
adiabatic couplings. These couplings are strongest near 
avoided crossings between adiabatic electronic potentials 
and lead to inelastic losses. For example, a transition 
from the M = 4 to M = 3 curves in Fig. |l| results in de- 
polarization and conversion of internal energy to kinetic 
energy. 

First we will assume that adiabaticity holds. We relax 
this assumption below. The zero-energy s-wave (I = 0) 
scattering length of polarized metastable 3 P2 m = 2 
atoms is then given by the scattering length of the Uq(R) 
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FIG. 4: Scattering length of the Uo(R) potential of spin po- 
larized Z P2 atoms as a function of magnetic field. Values for 
24 Mg, 40 Ca, and 88 Sr are shown. The arrows show the mag- 
netic fields for which zero-energy resonances occur. 



component of U2.2 (R, 9b)- This length for three alkaline- 
earth species as a function of magnetic field is shown 
in Fig. U For B < 40 mT the scattering lengths are 
large and positive and a Bose condensate would be sta- 
ble against collapse. Uncertainties in QQ and DD coeffi- 
cients do not change the picture significantly. 

For B > 70 mT the scattering lengths exhibit singu- 
larities, where a is infinite. These resonances are due the 
appearance of the first bound state in Uq(R). The inset of 
Fig. H shows that for increasing B the well becomes wider 
and deeper, supporting an increasingly larger number of 
bound states and thus extra resonances. 

Magnetic field control of the sign and magnitude of 
the scattering length near the resonances appears likely 
for these systems. Magnetic field-induced resonance be- 
havior in scattering lengths has been observed for alkali- 
metal atoms. There, however, the singularity is a conse- 
quence of Feshbach resonances rather than through the 
appearance of extra bound states in an adiabatic poten- 
tial. 

Inelastic processes — The scattering lengths have been 
determined assuming that nonadiabatic couplings and 
thus losses are negligible. We verify the applicability 
of the adiabatic approximation by estimating loss rates 
from a two-channel curve-crossing model. The model 
has an attractive van-der-Waals potential in an incom- 
ing channel, \i), of two spin polarized atoms and I = 
that crosses a repulsive C5/R 5 exit channel, |e), of two 
atoms with 777,1+7772 < 4 and I = 4. The coupling, C' 5 /R 5 , 
is due to the QQ interaction. More precisely, 
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where [i is the reduced mass of the dimer, A oc B is the 
difference in dissociation energy of the two states, the 
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exit state has a I — 4 centrifugal potential, and Cg is 
defined in Eq. [|. The two states cross at near Rb- 

The i?-dependent eigenvalues of V^ch are our model 
adiabatic potentials. In fact, the top-most adiabat plays 
the role of the L = component of the U22 adiabat. It 
has an attractive van der Waals tail and a repulsive wall 
that is to be compared to the "Zeeman-van der Waals" 
well of Fig. |. 

Loss can occur when the atoms reach short-range R on 
the incoming channel or when the atoms are reflected into 
the repulsive exit channel. We numerically solved the 
two coupled Schrodinger equations assuming a hard-wall 
placed inside Rb and parameter values valid for stron- 
tium. This approach has no losses due to short-range 
inelastic processes. Nevertheless, we feel that the model 
gives an order of magnitude estimate of the loss rate. 
The calculations show that for B from 1 mT to 5 mT 
the loss rate K into the exit channel rapidly rises from 
~ 10~ 13 cm 3 /s to 10~ n cm 3 /s due to Wigner-threshold 
effects in the exit channel. For large fields the rate de- 
creases slowly and exponentially with B and for B 500 
mT has dropped to 10 -12 cm 3 /s. 

The loss rate and the scattering length show a reso- 
nance similar to the resonances shown in Fig. |] when an 
extra bound state appears in the upper adiabat. For col- 
lision energies below 1 /iK the loss rate approaches zero 
in a narrow w 1 mT region above the resonance location 
where the scattering length is large and positive. For col- 
lision energies above 100 /iK the resonance is broadened 
and will be hard to observe. The position of the reso- 
nance and the corresponding loss rate will change when 
all coupled channels are included. However, the reso- 
nance and a (partial) drop in the loss rate is expected to 
survive. 



Consequently, we have identified magnetic field regions 
where losses are smaller than or on the order of 10~ 13 
cm 3 /s, i.e. for B < 1 mT and for magnetic field strengths 
just above the resonance. For several Tesla fields the 
loss rate will be small as well. At a rate constant of 
10~ 13 cm 3 /s the imaginary part of the scattering length 
hK/ (2/f)p0|, is orders of magnitude smaller than the real 
part of the scattering length shown in Fig. |ij The colli- 
sion is adiabatic for these field regions and the scattering 
length is meaningful. Moreover, these rates lead to sam- 
ple lifetimes of 1 s at densities of 10 13 cm~ 3 . 

Conclusion — We have calculated the interaction po- 
tentials of metastable alkaline-earth dimers in the pres- 
ence of a magnetic field. A purely long-range interaction 
potential is found to determine the scattering length be- 
tween two spin polarized 3 J*2 alkaline-earth atoms. A 
magnetic-field-induced resonance is observed, where the 
interaction potential supports an extra bound state. Ref- 
erence [|Io| predicts similar long-range states for polar 
molecules in an electric field. Existence of such molecu- 
lar states should be a general property of colliding atoms 
or molecules with an anisotropic interaction potential in 
an external field that split a degeneracy. 

Although full-scale multi-channel calculations of mag- 
netic trap losses and scattering lengths are desirable, a 
two channel model indicates that the inelastic loss rate of 
the collision between two spin polarized atoms are small 
for small magnetic field strengths and above the field- 
induced resonance, where the scattering length is posi- 
tive. We uncovered several unique features of collisions 
between metastable alkaline-earths, which may offer new 
insights into the physics of ultracold quantum gasses. 

This work was supported in part by the National Sci- 
ence Foundation and the office of Naval Research. 
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